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Abstract 



We investigate the effects of thermal interactions on tracking models of quintessence. 
We show that even Planck-suppressed interactions between matter and the quintessence 
field can alter its evolution qualitatively. The dark energy equation of state is in 
many cases strongly affected by matter couplings. We obtain a bound on the coupling 
between quintessence and relativistic relic particles such as the photon or neutrino. 



*hsu@duende. uoregon.edu 
Tbmurrayl@dar kwing.uoregon.edu 



1 



1 Introduction 



Recent evidence [1] suggests that a large fraction of the energy density of the universe has 
negative pressure, or equation of state with w = p/ p < 0. One candidate source of this 
dark energy is a slowly varying and spatially homogeneous scalar field called quintessence 
[2]. Because the dark energy redshifts more slowly than ordinary matter or radiation, it 
appears that the ratio of energy density of quintessence to that in ordinary particles must 
be fine tuned to a specific infinitesimal value in the early universe in order to explain its 
current observed value. One class of models that ameliorate this problem describe tracker 
fields [3] whose evolution is largely insensitive to initial conditions and at late times begin 
to dominate the energy density of the universe with a negative equation of state. Tracker 
models have difficulty producing consistent with observational data: they generally imply 
wf > -0.7 , whereas WMAP implies wf < -0.78 (95% CL) [1]. (Here, effective means as 
measured observationally, so integrating over redshifts less than of order 10 3 .) Nevertheless, 
they provide an interesting class of models describing dark energy as a slowly evolving scalar 
field. An alternative class of models, which avoids the extremely flat potentials required at 
late times of tracker models, utilizes nonlinear field oscillations that exhibit w < [4]. 

Tracker models generally require only a single adjustable parameter. Once this parameter 
is appropriately chosen, a wide range of initial values of the tracker field, 0, and its derivative, 
0, result in similar values of its energy density today. This is due to an attractor-like property 
of the tracker equations of motion. 

In this paper we investigate the effects of interactions between the quintessence field, 0, 
and ordinary matter particles in the early universe. It is important to note that while the zero 
temperature potential may be fine tuned in order for the evolution of to have the attractive 
properties mentioned above, the same may not be done with the finite temperature effective 
potential. That is to say, once the form of the renormalized zero temperature potential 
is determined, no additional freedom remains to fine-tune away unwanted thermal effects. 
Therefore, such effects must be considered. We expect thermal interactions to be at least of 
gravitational strength (even in the case where is a "hidden sector" field). At minimum, 
quantum gravity is likely to produce interactions of the type [5] 
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■A , (i) 



where Mp is the Planck scale, and are terms in the standard model lagrangian, including 
for example 

F% , F^ , $W , • • • 
where F is the field strength of any gauge field (including the photon, but not excluding 
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gluons or the W or Z) and ip is any fermion field from neutrinos to the top quark. Even 
if were a pseudo-Goldstone boson [6], it would be surprising not to find at least Planck- 
suppressed violations of the resulting — > + constant symmetry. String theory, for 
example, is believed to not exhibit any exact global symmetries [7]. Previous constraints on 
certain $ are quite strong, where the coupling is to the photon or gluon [6] . However, some 
fa could be much larger, such as when the interaction is with the W, Z or even a neutrino 1 . 

In the early universe matter particles are in thermal equilibrium, and the interactions in 
(1) produce a thermal mass for of the form 



where T is the temperature. If the thermal degree of freedom is massive, the thermal effect 



note that the thermal effects of interest here are in addition to any quantum corrections to 
the effective potential resulting from the interactions between matter and quintessence (see, 
for example, [9]). In general, the bare parameters of quintessence models must be fine-tuned 
in order to obtain potentials of the necessary type. We assume here that this fine-tuning is 
achieved (whatever its consequences for the plausibility of the model) and focus on thermal 
effects which must also arise. 

Although may not itself be in equilibrium, nevertheless its dynamical evolution will 
be affected by these thermal interactions, just as for the axion field near the QCD phase 
transition [10]. We can derive the correction (2) to the effective potential for as follows. 
Let the cold field be a static, external source for a Euclidean path integral describing the 
thermal degrees of freedom. The timelike boundary conditions for the path integral have 
period given by the inverse temperature. Performing the integration over the thermal fields 
yields a contribution to the effective potential for 0, and the usual perturbative analysis 
identifies the leading effect to be the thermal mass term in (2). In this calculation, we 
need never assume that itself is in thermal equilibrium, yet its effective potential receives 
temperature-dependent contributions. 

In what follows we will examine how the contribution of (2) to the tracker potential 
modifies its evolution. We can make a simple argument for why (2) is non-negligible at late 
times. At late times, the quintessence field must have a very small mass: l / "(0) 1 / 2 <Hq^ 
10 -33 eV, and contribute of order closure density to Q: V(<f>) ~ (10 _3 eV) 4 , which implies 
that ~ M P . This means that the mass term in (2) can be roughly the same size as V(4>), 
up to powers of 

1 Direct coupling to relic neutrinos has been considered previously [8]. 
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goes to zero exponentially as e 



l T when the temperature drops below the mass m. We 
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At early times, (2) also affects the evolution in many cases. Suppose the tracker potential 
is given by V(<f)) = M /+4 _/ , where / > 4. Then V(<f)*) and (2) are comparable at the 
minimum of the combined potential: 

where the potential energy density is roughly 
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In many cases, oscillates about the temperature-dependent minimum 0*. The oscillation 
energy redshifts faster than the potential energy at the minimum, V* ~T 4 '/('+ 2 ), so simply 
tracks 0* with oscillations that decrease in amplitude over time. Interestingly, K redshifts 
exactly as the energy density of the tracker solution [11] (assuming radiation domination; 
during a matter dominated epoch V* redshifts somewhat faster than the usual tracker energy 
density). This means that thermal effects will keep and its energy density near their 
desired values, even though the physics responsible is very different. When the thermal term 
eventually either disappears due to the crossing of a particle mass threshold, or becomes 
negligible due to redshift, will merge back to a tracker solution. 



2 Evolution results 

We assume a spatially flat Robertson- Walker universe, with metric ds 2 = dt 2 — a 2 (t)dx 2 . 
The evolution of a scalar field minimally coupled to gravity in this spacetime is given by the 
Klein-Gordon equation: 

+ 3#0 + V'((f)) = , (5) 

where a dot denotes the derivative with respect to cosmic time, and a prime denotes the 
derivative with respect to 0. The evolution of the scale factor is governed by the Friedmann 
equation: 

where if z denotes the redshift, then p m = p c Q m (l + z) 3 , p r = p c fi r (l + z) 4 , and p<f, = 
^0 2 + V(4>). Here the subscript m refers to both baryons and cold dark matter, and the 
subscript r refers to both photons and neutrinos. If the universe is spatially flat, then it will 
always be the case that Q m + Q r + = 1. Observational data [1] currently favor fl m ~ 0.3, 
tt r ~ 1(T 4 and ~ 0.7. 
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Equations (5) and (6) were integrated numerically for a wide range of 0j and 0j from an 
initial redshift of Zi = 10 28 (temperature ~ 10 16 GeV), which might plausibly correspond to 
the end of inflation. Motivated by the arguments of the previous section, we took V((b) to 
be 

V(cf>) = M 4+ <0~' + {^j 2 T 4 , (7) 

where / > 4, (3 is a free parameter, and M is constrained such that ~ 0.7. For example, 
for I = 6 and f3 = 0, M ~ 4.7 x 10 6 GeV. In this simulation, T(z) = tfp~ r = ^A(l + z) 
(we are not precise about the number of relativistic degrees of freedom). It is worth noting 
that, since the second term in (7) arises due to a loop effect, it should actually appear in 
(7) multiplied by a constant of order 1CT 1 . In the absence of this factor, the quantity (3 in 
(7) differs somewhat from the $ in (1). This consideration, however, has no effect on the 
qualitative picture described below. 

Based on our simulations, we make the following observations. 

For (3 = 0, tracking occurs for a large range of initial conditions in and 0, as described in 
[3]. In particular, for / = 6, if starts from rest, any 0j in the range 10" 18 Mp < 0« < 10~ 2 Mp 
will be on track by today. In general, the limits for 0« that will be on track by today, assuming 
<pi — 0, are found by solving p^ = M l+4 (f>~ 1 for 0j, where is the initial energy density in 
0. By noting that M ~ (p < p M l P ) 1 ^ 4+l \ where p^, is the present energy density in 0, it is 
straightforward to see that these limits on 0j depend on /: 

t.^ Mp (&A l . (8) 

The minimum value of 4>i that will be on track by today, 4>^ m i ni is then simply found by setting 
p^i equal to its maximum value. For an initial redshift of 10 28 , this corresponds roughly to 
Pfa ~ 10~ 4 pBij where ps% is the energy density of dominant background component, radiation 
at this redshift (p r {z = 10 28 ) ~ 10 61 GeV 4 ). Similarly, the maximum value of <pi that will 
be on track by today, <pi^ maX) is found by setting p^j equal to its minimum value, which is 
roughly the background energy density at equality, p eq ~ 1CT 37 GeV 4 . 

For f3 7^ 0, (pi^min is essentially unchanged because the first term in (7) is dominant for 
P0 i ~ 10~ 4 p_Bi, unless (3 is made very large (f3 ~ 10 20 ). (3 this large will not be discussed 
further in this paper. For sufficiently small (3, (pi tmax is also left unchanged. Let (3 C = 
(10- 4 ) 1 /2(p e(? /p^)V „ io(ii-20A Then for (3 = (3 C and 0, = 0^, P(j>i = l(T 4 p K . But 
P4>i < lO^pBi, and p oc f3 2 (b 2 for > 0*. Therefore, for (3 > f3 C) (b^ max oc l/f3. The net 
result is that the range of 0j that will be on track by today (with 0j = 0) is independent of 
(3 for (3 < (3 C and goes like 1/(3 for (3 > (3 C . 

In addition to affecting the range of 0j that track, the choice of (3 qualitatively affects 
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the dynamics of 0. Note that 10 -2 < f3 c < 1. In examining the dynamics of (b there again 
seems to be a critical value of (3. Although this critical value seems to be ~ 10 _1 , it is not 
clear whether it is equal to f3 c . As mentioned above, a factor of order 10 _1 was not included 
in (7). Therefore, in what follows reference will be made to /3 C , which is meant to indicate a 
(3 in the range 10 -2 < (3 < 1. 
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FIG. 1. The evolution of p^ is shown for (3 = 3, I = 6 and an intermediate p^ (<pi <C 0*j). p^ for 
4>i = <f>*i (corresponding to the tracker solution for (3 = 3) is also plotted for reference, as are p r 
(medium gray) and p m (light gray). Note that p^ cannot decrease below the the energy density of 
the tracker solution (it is always larger than V*) and subsequently freeze as it does for (3 = 0. Also, 
note that, at high redshift, small oscillations can be seen in p<p corresponding to oscillations of (f> 
about 4>* (as in Figure 2). 




FIG. 2. The late time behavior of the tracker solution is shown for (3 = (lower curve) and (3 = 10 2 
(upper curve) with I = 6. Again p r (medium gray) and p m (light gray) are plotted for reference. 
Note that while p${(5 = 0) has a rather shallow slope today, p</,(/3 = 10 2 ) cannot because p^ > V*, 
and V*((3 = 10 2 ) > p m . Also, note that the redshift at which p^ begins to dominate depends on (3. 

For (3 < fie, the behavior of 4>(z) is essentially just that described in [3], with the one 
additional constraint that p<f,(z) > V*(z), at all times (see Figures 1 and 2). For (3 > (3 C) 
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<j)(z) oscillates about (f>*(z) (see Figure 3). The period of these oscillations in <f)(z) decreases 
exponentially with the scale factor, while the amplitude decreases monotonically, but not 
exponentially, per se. As either side of the potential (7) are made steeper; i.e., either I or 
f3 are increased, the period of oscillation decreases, as does the amplitude. The fact that 
these oscillations have been damped out by today and that V* redshifts at the same rate (for 
RD) or faster (for MD) than the tracker solution is what allows for tracking to occur, even 
with (3 > f3 c . Adding extra terms to (7), which cause V* to redshift slower than the tracker 
solution during RD will be discussed below in Section 3. 
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FIG. 3. The evolution of <fi is shown for f3 = 3, I = 6 and an intermediate (<fii <C 0«). 4>* is 
also plotted for reference. Note that <fi oscillates about 0* at large redshift; follows 0* closely for 
most of its evolution, but for small redshift it begins to fall behind. 

For (3 < /3 C , the scale M in (7) is simply equal to its value for f3 — 0; i.e., M ~ 
{p < f )0 Mp) 1 ^ 4 ' +l \ For (3 > [3 C) M slowly decreases as (3 is increased, e.g. for I = 6 and (3 = 0, 
M ~ 4.7 x 10 6 GeV, whereas for I = 6 and (3 = 10 3 , M ~ 1.0 x 10 6 GeV. 

We define the scalar field equation of state: 

_ 1 + z dp^ 

The definition (9) coincides with the usual definition of the scalar field equation of state: 

as a consequence of energy conservation. However, because we did not take into account 
the back reaction on matter and radiation of the interaction, there are regimes in our 
simulation where (9) differs from (10). To be precise, we have treated ordinary matter as a 
thermal background, and have not accounted for energy flowing from into the heat bath. 
This is generally a negligible effect, except when the dark energy density is large and 
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is of order Mp, which can occur at late times. The evolution of can strongly influence 
the thermal matter (for example, changing the coefficent of its kinetic term), which we 
have not accounted for. The late time behavior of our simulations at large f3 is therefore 
only qualitatively and not quantitatively correct. Observational probes such as WMAP are 
sensitive to the way the dark energy density redshifts, and hence constrain (9). 



W0 




FIG. 4. The evolution of is shown for (3 = 10 2 , / = 6 and an intermediate p^i (0j <C 0«). Note 
that the (very rapid) oscillations in uu (corresponding to oscillations of about 0*) are completely 
damped out very early in the evolution of 0, and that ~ for most of the evolution of 0. Also, 
note that today ~ —0.13, compared to the (3 = value of ~ —0.4. 
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FIG. 5. The present value of w^ is plotted versus log(/3) for I = 6. Note that for (3 > 10, w^ Q 
begins to deviate from its (3 = value. For (3 > 10 3 , w^o ~ 0. 

F° r P ^ Pc there are oscillations in w^, corresponding to the oscillations seen in (see 
Figure 4). Again, the period of these oscillations decreases as ft is increased, and the period 
decreases exponentially with the scale factor. By increasing /, the present value of is 
driven toward zero (from below). By increasing (3 (beyond ~ (3 C ), w<f, is also driven towards 
zero (see Figure 5), as the evolution of becomes controlled by 0*, and the energy density 
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by V*. In the extreme limit tracks 0* closely and its energy is almost entirely potential, 
rather than kinetic. Note that — 1 < < at late times for all / and (3. 

For (3 — and / = 6, w<f, ~ —0.4. In general, if a sum of inverse powers of are allowed, 
then w'f > —0.7 [3]. This is the effective equation of state measured by supernovae and 
microwave background experiments, which integrate over a (potentially) varying w</,. This 
bound represents the best case scenario for (3 = tracker models, in that it is most consistent 
with observational data for u^ ff , which put u^ ff < —0.78 (95% CL) [1]. The inclusion of the 
second term in (7) with sufficiently large (3 {(3 J> (3 C ) results in ~ 0* throughout most of 
its evolution. For I = 6, this yields an equation of state, w ~ (since V* ~ T 3 ), that is even 
further from the observational bound (see Figure 5). At very late times, does not increase 
sufficiently rapidly to stay near 0* (see Figure 3). Instead, it rejoins a tracker solution and 
its equation of state reverts to one in which the energy density redshifts more slowly than 
V* (see Figure 2). 

3 Discussion 

Our analysis shows that the evolution of the tracker field depends quite sensitively on its 
interaction with thermal matter, even when the strength of the couplings is as small as 
one would imagine they may possibly be; i.e., Planck-suppressed. When (3 is larger than 
of order unity, there is a tendency for the evolution to be controlled by that of 0*, once 
initial oscillations have damped away. By a lucky coincidence, the redshift of V* ~ T 4l ^ l+2 ^ 
is the same as that of the tracker solution (during radiation domination), so that can 
rejoin a tracker solution at late times. The most dangerous possibility (which is realized 
when (3 > (3 C ) is that is still following 0* at late times, in which case its equation of state 
will be far from the observationally favored = —1. For general I, the equation of state 
obeyed by K is = (I — 6)/3(Z + 2), which is never consistent with observational bounds 
for I > 4. We have checked that the behavior described above is qualitatively similar when 
higher order terms such as 



are included in the potential. 

Hence, we conclude that there are stringent limits on the coupling between the tracker 
field and any particles which are still relativistic today, such as the photon or neutrinos. Such 
limits cannot be avoided through fine tuning of the finite temperature effective potential; 
once the zero temperature potential has been computed the finite temperature effects are 
determined. Interactions which are more than roughly two orders of magnitude stronger 
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than Planck- suppressed lead to a problematic equation of state. Couplings of the tracker 
to heavy particles, which freeze out at T ~ m, may alter the tracker evolution at early 
times, but do not affect the observed dark energy equation of state and are hence poorly 
constrained. The best hope of directly detecting the quintessence field may be through its 
interaction with massive particles. 

Finally, although our analysis has focused on tracker models, similar results apply for 
any quintessence model in which the field is today slowly evolving in a very flat potential. 
As we argued in the introduction, in any such model the value of must be of order Mp 
today, which means that thermal terms such as (2) or (11) will be important for sufficiently 
large (5. Large couplings to relic particles such as neutrinos can be ruled out as they lead to 
a problematic equation of state. 
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